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HYBRID CALIBRATION OF TISSUE DENSITIES IN COMPUTERIZED 

TOMOGRAPHY 

Background of the Invention 

Field of the Invention 

[0001] The present invention is in the field of medical imaging using computerized 
tomography (CT), and in particular, is directed to a method to determine tissue densities 
in the body of a subject and to provide calibrated display of images. 

Description of the Related Art 

[0002] CT scanners have become a major diagnostic modality in modern medicine 
and are widely used for many types of exams. Most exams are performed by subjective 
viewing of the cross-sectional images on either film or electronic displays. This 
subjective viewing makes use of assumed quantitative image pixels, which define 
boundaries of tissues, organs or foreign masses, and subjective discrimination of tissue 
types by density differences. Identification of diagnostic details is fundamentally 
dependent upon the detection of image detail edges* 

[0003] Measurements of true tissue densities in the living subject have many 
diagnostic benefits, in addition, to bone densitometry. Several new and promising 
measurements include lung nodule density, cardiac calcifications, aortic calcifications, 
soft plaque, fat measurements, BMI, lung volume and density, liver iron content, and the 
like. Knowledge of true tissue densities will allow diagnostic analysis of images not 
currently possible. Absolute change in CT numbers may allow new diagnostic criteria. 
Emphysema, tissue fat content, calcifications, liver iron build up, and the like could be 
determined from the calibrated data, thus, adding a new dimension to CT interpretation. 
[0004] Radiologists routinely make subjective, and even quantitative measurements 
of foreign masses, tissues or organs by manually placing cursors to define the 2-D extent 
of the target. If the window and/or level (brightness and contrast) are changed in the 
display, the apparent size of the target changes because the boundary is not discrete and is 
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moved in or out of the display range. The measured object size is, thus, frequently 
inaccurate, and will vary from operator to operator and from scanner to scanner 
depending on the display conditions and scanner properties. 

[0005] CT images are filmed by a technologist or other operator and are recorded on 
standard x-ray film for light box viewing. The size, and apparent density of target 
objects, and foreign masses depend on the window and level settings. The window/level 
settings are subjectively set for filming and display of a particular image. In addition, the 
process to set and adjust the window and level requires operator time and is currently 
very inefficient. Electronic image data are frequently erased, and only the films retained 
for the medical records. Later viewing is limited to the subjective display and/or the 
filming levels previously chosen by the operator. 

[0006] The foregoing discussion is based on the assumption that pixels and/or voxels 
of the image are a representation of the true underlying density of the target tissue. 
Although this assumption is roughly maintained due to the scanner being calibrated to 
water and air, it is sufficiently inaccurate that many quantitative measurements cannot be 
made with even the best modern scanners. 

[0007] There has been significant, recent interest in quantifying coronary calcium, as 
well as calcifications in the aorta, lungs, breast, and carotids. It is desirable to provide 
improved calibration methods for all the tissues of the body. Whole body CT scanning is 
growing rapidly in use. The entire torso is scanned creating many thin slices for analysis 
and viewing. The radiologist attempts to subjectively analyze many structures from 
many images, which is laborious and very time consuming. Measurements of densities 
and volumes of many organs are of interest, including heart, lung, liver, kidneys, prostate, 
thyroid, pancreas, and the like. Quantitative measurements of calcifications, of areas and 
volumes, as well as standardized viewing and filming of images, all require improved 
calibration methods. 

[0008] CT scanners have been used as quantitative instruments for bone density 
measurements in quantitative computerized tomography (QCT) by the use of calibration 
phantoms. More recently, fast CT scanners, such as the Imatron EBCT and GE Light 



Speed, have been used for coronary calcium analysis with or without phantom 
calibration. Several calibration approaches have been used in QCT bone densitometry 
including simultaneous phantom calibration with bone and tissue equivalent phantoms, 
non-simultaneous calibration with more anthropomorphic phantoms, non-phantom 
calibration using histogram analysis of fat and muscle regions, simultaneous phantom 
calibration with blood sample corrections, and dual energy calibration to correct for 
vertebral fat in bone density measurements. These approaches have been developed 
specifically for and used for QCT bone densitometry of vertebral trabecular bone. 
[0009] CT numbers, (Hounsfield Units, HU), are estimates of the attenuation 
coefficients of tissue relative to water as the calibration reference material. However, CT 
numbers fail to be truly quantitative for several reasons. For example, the tissue 
attenuation coefficients are photon energy dependent, and the x-ray beam energy spectra 
are not measured or known for individual patients. Further, there exists many beam 
energy spectra in each CT slice, i.e*, a unique spectrum for each path length through the 
patient, and seen at a particular detector element, and a unique spectrum for each view 
through the patient. The beam spectrum changes with the thickness and composition of 
tissues in the path length. The quantities of fat, soft tissue, air, and bone vary with each 
projection. X-ray tube filtration to shape the beam intensity also changes the beam 
spectrum resulting in variation in CT numbers based on locations within the field of 
view. Image processing software and current beam hardening corrections have as an 
objective to improve subjective image quality, and do so, often, at the expense of 
quantitative information. CT number calibrations and beam hardening corrections are 
based on idealized phantoms, which are often circular in shape and composed of water, 
plastics, or other synthetic materials. These differ significantly from the shape and 
composition of real patients. CT numbers at the edge of the field of view, where a 
calibration phantom would be placed, are different from those inside the patient. This 
produces errors in calibration since the phantom can never be placed inside the body 
cavity. CT numbers vary through the image on each slice, and are dependent on table 



height, position in the beam, slice thickness, field of view, and sometimes even the time 
of day as the scanner warms up. 

[0010] Many diagnoses are based subjectively on perceived tissue densities and 
regional changes in density as demonstrated by the CT numbers in the image. Results 
currently are independent of patient variability and CT equipment* Standardization and 
calibration of the CT numbers across different patients and CT scanners will aid in 
interpretation of many conditions. 

[0011] It is frequently desirable to make quantitative measurements from both two- 
dimensional (2D) and three-dimensional (3D) data sets in medical imaging. Accurate 
measurements of organ or tumor volumes and cross-sectional areas of various biological 
details, such as blood vessels, all have potential medical diagnostic value. Quantification 
of vascular calcium and micro calcifications throughout the body is valuable in 
cardiovascular disease and breast cancer detection, for example. All of these tasks use 
gray scale, voxel-based data. The identification of the edge of a target region may use 
any of several edge detection algorithms, such as the Sobel operator in either 2D or 3D 
space. This prior art has used image voxels as outputted from the imaging device, i.e., 
CT scanner, digital radiography apparatus, magnetic resonant scanner (MRI) or 
mammography system. In all cases, the image data was not calibrated. Since the image 
gray scale values vary with the imaging conditions and subject properties, the definition 
of an edge also varies. 

[0012] In some cases, the diagnostic detail is defined by a pre-selected threshold 
value, i.e., if the target element equals or exceeds the threshold value, the detail is 
counted as a positive diagnostic find. Coronary artery calcifications are a notable 
example. With currently available CT scanners, calcifications that exceed either 130HU 
or 90HU are counted as positive finds. The Hounsfield units (HUs) are known to vary 
with scanner type, x-ray beam energy, reconstruction software, patient size and 
composition, and the like. As a result, the threshold value varies depending on these 
conditions. A positive calcification find is thus different for a small female versus a large 
male. If a patient is scanned on one scanner and scanned later on a second scanner for a 



follow-up exam, the results will be different. The diagnostic results are therefore 
dependent upon several variables of the imaging systems, as well as being dependent on 
the patients. 

[0013] The use of external calibration phantoms containing bone equivalent samples 
have been used for some time in QCT Bone Densitometry. Such phantoms have greatly 
aided the standardization accuracy and reproducibility of bone density measurements. In 
this case, however, the target tissue, bone, is large, of a high density much larger than soft 
tissue, and located relatively close to the calibration phantom. 

[0014] The use of external calibration phantoms has only recently been attempted 
with coronary calcium quantification. Calibration phantoms have not been used for soft 
tissue density measurements or for physical dimensional measurements. One of the 
problems which exists with quantitative CT relates to the variation of image gray scale 
values throughout the area of the image. The same tissue type located in one location 
within the body may produce a different Hounsfield unit value versus a different location. 
That is, the image is not homogeneous throughout. Not only does the image vary in 
intensity, but it also varies in effective beam energy. As a result, no one unique 
calibration curve is available for each CT slice or for a complete digital 2D radiograph. 
The situation is complicated by being dependent upon the position of the object within 
the scan field, device servicing and calibrations, and x-ray tube wear. 

Summary of the Invention 

[0015] Unlike prior art devices in which image data was not calibrated, the disclosed 
invention presents a method to calibrate the image data such that a calibrated image detail 
edge can be quantitatively defined, located and relocated on follow up images. The 
reproducibility of physical measurements of target volumes, areas, and distances is 
greatly improved. 

[0016] The embodiments of the present invention provide a method to calibrate the 
image to produce consistent and standardized results, independent of the measurement 
conditions discussed above (e.g., scanner type, beam energy, reconstruction software, 
patient size and composition, and the like). Clinical results will be more consistent, more 



accurate, and more precise. The ability to follow patients on follow-up exams to monitor 
disease progress will be greatly benefited. 

[0017] The density differences in Hounsfield unit values between the soft tissues of 
the body are small, compared to bone density, and they are more difficult to measure. 
Microcalcifications in the vascular system of the body produce small density differences 
as well, and are difficult to measure reproducibly. The embodiments of the present 
invention overcome the small differences in the density of soft tissue or 
microcalcifications to facilitate measurements throughout the body, both in areas near the 
external phantom and in areas distant from the external phantom to overcome the 
inadequacies of the conventional methodology of QCT calibration with an external 
phantom. In particular, the embodiments of the present invention enable a method to 
quantify x-ray images that uses a hybrid calibration method that overcomes these 
limitations. 

[0018] Filming and viewing CT images can be automated and presented in a 
quantitative way, which further aids the interpretation. Standardization of filming 
windows and levels to an absolute scale based on true density provides consistent 
viewing conditions. Accurate and reproducible filming could be automated and 
standardized to quantitative values consistent through the medical community. 
[0019] The body is always in a state of homeostasis, attempting to keep the body in a 
static and stable condition. Indeed, if some tissues of the body change by even small 
quantities, life may be challenged. Blood is the fluid of life, and must be maintained 
within relatively narrow ranges. It is therefore a very consistent substance, both in 
individual patients, and from patient to patient. The literature on blood composition from 
the standpoint of effects on physical density shows that blood density and the resulting 
CT HU values do vary with the hematocrit, red cell volume, iron content, blood fats, and 
hydration, but these variations are small compared to the required calibration accuracy. 
This analysis indicates that these variations are acceptable for the disclosed invention. 
Biologic variations in blood density are expected to generally produce a range of 
variation of ±2 HU and to rarely exceed ±4 HU. 



[0020] Since blood is always homogeneously mixed and consistent throughout each 
individual body, the blood provides for a unique and consistent in vivo calibration tissue. 
Since blood in the great vessels and heart is centrally located within the body, the 
position of such blood is ideal as an internal surrogate calibration tissue. The 
embodiments in accordance with the present invention provide a method that uses blood 
for tissue calibration in x-ray imaging systems, both when used with an external 
calibration phantom and when used without an external calibration phantom. 
[0021] Air calibration is frequently used in CT scanner calibration along with water. 
The Hounsfield unit value of air is defined as the minimum CT density, usually -1000 
HU. Air density is used for quality assurance and routine calibration-file setup of the 
scanner. Air calibration has not been used as a calibration reference with individual 
patient scans. It can be assumed that air is a consistent and reproducible substance. 
Scatter radiation degrades images by adding a DC image component, which reduces 
image contrast, signal-to-noise ratio, and dynamic range. The measurement of tissue 
densities and target edges are degraded. The air present within the esophagus provides an 
internal air reference for calibration in cardiac and chest imaging. Bowel gas provides a 
potential air reference for the abdomen. The air adjacent, but outside the body in the 
environment, provides an external air reference. The embodiments in accordance with 
the present invention provide a method to use both internal air and external air as 
calibration references in tissue density measurements. 

[0022] Prior art has disclosed methods to use muscle and visceral fat with histogram 
analysis to provide calibrations without an external phantom for QCT bone densitometry 
(See, for example, U.S. Patent No. 5,068,788 to Goodenough). Muscle and visceral fat 
always have mixtures of both tissue types, muscle and fat, but the technique proposes that 
the histogram peaks can be made to represent each tissue type as calibration points. It 
would be helpful to have a tissue which contains fat without muscle intermixed. This 
does occur, but, unfortunately on the exterior circumference of the body as subcutaneous 
fat. 



[0023] Segmentation techniques are now available to isolate these fat tissues. The 
thin or not so thin layer of subcutaneous fat can be defined and isolated with modern 
segmentation techniques. The location of the fat tissue, within the scan field and relative 
to the target organ, can be measured. Since the fat is distributed around the subject, the 
fat can be used to make corrections for beam hardening and scatter. The present 
invention discloses a technique for using subcutaneous fat as a known reference tissue for 
calibration in tissue density measurements with a CT scanner. 

[0024] One aspect of the embodiments in accordance with the present invention is a 
method to calibrate CT images. The method utilizes both an exterior fixed calibration 
phantom with known samples of tissue equivalent materials, and which also uses an 
interior sample of tissue of the patient. The exterior calibration phantom has multiple 
samples of varying densities, varying compositions, or varying densities and varying 
compositions, which provide a common standard of reference that is consistent for all 
scanners. The phantom sample readouts provide a regression equation which has a slope 
that provides a measure of the effective beam energy of the CT scanner for that particular 
patient. The phantom samples run substantially along the torso of the patient, and the 
samples are included in each CT slice with simultaneous scanning. The calibration 
offset, the y-axis intercept, is determined from the patient's own tissue, preferably heart 
tissue and/or blood. The total heart is segmented in 3D space, and a best representation 
of the average voxel HU value is determined. Calcifications and fat are removed from 
the volume by thresholding and by histogram analysis. The calibration makes the 
assumption that blood density is consistent, and the same for all people. Although blood 
density does vary with hematocrit, blood cell volume, and iron content, these variations 
are relatively small and acceptable. Human blood, and heart tissue are remarkably 
similar in all subjects. In addition, the embodiments of the present invention make use of 
the assumption that blood is homogeneous and has a constant density in every subject, 
which allows calibrations using the heart and great vessels. 

[0025] A hybrid calibration equation is created for each CT slice, and thus each 
region of the body is calibrated with data from that region. The blood and heart tissue 
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become an additional in vivo phantom sample, which is homogeneous and is dispensed 
throughout the body. The hybrid calibration is applied to all image pixels of all slices to 
generate a new data set, which is now calibrated and standardized. 
[0026] Absolute thresholds and standard edge detection algorithms are applied to the 
calibrated images to measure volumes, areas, and distances. Organ or tumor mass can be 
readily measured in a reproducible manner, for example: Absolute and quantitative 
display ranges, window and level, are set as default values on electronic displays and for 
filming. Optimum and standardized viewing conditions can be readily maintained for 
more consistent diagnosis. 

[0027] The calibration methods disclosed herein can be readily applied on all CT 
scanners and on all tissues of the body with accurate and reproducible results while 
reducing operator time. Any tissue in the body can be calibrated including lung, heart, 
kidney, bone, calcifications, plaque, and the like. 

Brief Description of the Drawings 

[0028] Figure 1 illustrates a cross-sectional image through the heart with an exterior 
calibration phantom, in accordance with an embodiment of the present invention. 

[0029] Figure 2 illustrates hybrid calibration curves at different body locations. 

[0030] Figure 3 illustrates a histogram of the HU values of the segmented heart 
volume showing the mode which is used for the calibration method of the present 
invention. 

[0031] Figures 4 illustrates the gray scale values at the edges of a tissue detail and the 
influence of display window and level on the location of the perceived edges and size of 
the detail. 

[0032] Fig. 5 shows histogram plots of 4 ROIs positioned in air outside and around 
the body of a patient. 

[0033] Figures 6A and 6B illustrate a flow chart that represents the computational 
procedure of the hybrid calibration method. 



[0034] Figure 7 illustrates a cross sectional image of the abdomen with reference 
calibration phantom in place. 

Detailed Description of the Preferred Embodiment 

[0035] A preferred embodiment of the present invention is described below in 
connection with the flow chart of Figures 6A and 6B in view of the illustrations in 
Figures 1-5. 

[0036] As illustrated by a step 110 in Figure 6A, a CT scan of the subject is first 
taken with the subject lying on the reference calibration phantom. Several CT images are 
taken in a short time period. Preferably, CT images are taken using a multi-slice spiral 
scan CT scanner or using the EBCT Imatron scanner, although any CT scanner can be 
used for many tissues. Scans of the heart for coronary calcium analysis will require fast 
scan times to stop the motion of the beating heart. 

[0037] After reconstruction of the CT images in the scanner computer, a first cross- 
sectional image of the subject is displayed in a step 112 in Figure 6A. Figure 1 illustrates 
a representative depiction of such an image taken through the heart showing coronary 
calcium 6. The heart tissue 2, the aorta 3, the lung 4, the chest wall and a reference 
phantom 5 are shown in Fig. 1. The reference calibration phantom 5 includes three 
representative samples (SI, S2, and S3) of varying densities. A prior art automated 
algorithm is used preferably to find the calibration phantom in the image as shown in 
Figure 1, and to place Regions-of-Interest (ROIs) within each phantom sample, as 
illustrated by a step 113 in Figure 6 A. An exemplary prior art algorithm is illustrated in 
U.S. Patent No. 4,922,915 to Arnold, which is incorporated by reference herein. 
[0038] The phantom samples contain known sample densities (Sd) and produce ROI 
mean CT numbers (HU S ) for each sample. For example, sample Si may preferably 
contain water density or 0 concentration of the calibration additive material. The 
calibration additive material is advantageously a calcium containing material such as 
calcium hydroxyapatite in a solid water equivalent synthetic material. 
[0039] In a step 114, the phantom measurements are expressed in a regression 
equation of the form y = I 0 + SxHU s . The slope, S, of the regression curves is related to 
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the effective beam energy used to make the image. For different KVps, filtration of the 
primary beam, and different patient sizes and compositions, the slope will vary. The 
intercept, I 0 , of the regression equation is a measure of the CT number of the water 
density sample. The slope, S, will vary with location in a given patient due to differing 
amounts and densities of tissue being present in the beam. Figure 2 is a depiction of three 
regression equations showing different slopes for 3 regions of the body. These regions 
could be different cross sectional areas and/or different body compositions of fat, muscle 
or bone. The regression slope will also vary with CT scanner model and calibration. 
[0040] The reference tissue region is next located in a step 115. The reference tissue 
may be blood, the heart filled with blood, subcutaneous fat, internal or external air. 
Preferably the heart and blood are used in coronary calcium analysis. Other tissue may 
be preferable in other exams. These tissue regions may be located by making use of their 
location and expected CT number range and their location relative to the location of the 
calibration phantom. The operator may place a cursor marker within the reference tissue 
as a seed point. The software then uses region growing techniques to find the region 
boundaries. 

[0041] After either automatic location or manual location, the measuring analysis 
steps can be applied. One or multiple reference tissues may be used in the calibration. If 
the reference tissue is the heart and blood, the algorithm identifies the blood filled heart 
which is surrounded by lung tissue. The location step is aided by first finding the aorta 
and the inferior margins of the heart, which may be positioned in contact or adjacent to 
the diaphragm, the liver or both. The operator may need in some cases to place a cursor 
marker to aid in this boundary distinction. 

[0042] In a step 1 16, the regression equation is applied to the entire image, or, in 
some cases, may be applied to only the reference tissue region plus additional 
surrounding tissue. This step provides an improved definition of the reference tissue 
boundaries. Alternatively, the regression equation may not be applied to the image or 
reference tissue region. The methodology will work without this optional step. 
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[0043] In a step 1 17, the regression equation is used to define a quantitative value as 
the specific edge. By using calibration, the boundaries of the reference tissue can be 
reproducibly located and relocated with fixed criteria. This edge definition and 
calibration is not to be confused with the calibration step which defines the edge of the 
target tissue that will be discussed below. 

[0044] The voxels within the boundaries of the reference tissue are next used to 
calculate a histogram in a step 118. A representative histogram of the heart as reference 
tissue is shown in Figure 3. 

[0045] The histogram is smoothed and a Gaussian curve is fit to the peak. Several 
neighboring slices may be scanned at this point in a step 119, The ends of the histogram 
may be filtered to remove unwanted voxel readings not representative of blood. Then, in 
a step 120, a background level due to scattered radiation or surrounding tissue, which can 
be distinguished from the reference tissue, may be subtracted. Fig. 5 shows changes in 
scattered radiation outside and adjacent to the body. ROIs were positioned in air and the 
histogram plotted of the pixel values of the ROI. With correct air calibration and no 
scattered radiation, the ROI should produce HU values around - 1000. The mode of the 
plots produce values of -964, -980, -987 and -994. These measurements show 
differences from expected values of 6 to 36 HU, which differences depend on location. 
Similar measurements can be made on air within the body. 

[0046] The mode of the Gaussian fit is next calculated in a step 121. The calculated 
mode is next combined with the phantom regression equation to create the final hybrid 
calibration equation in a step 122. The calibration equation is an expression of the form: 

T d = (HU b -HU' b ) + SxHU, 

where Ta is the tissue density of the target tissue or organ to be analyzed. Tb may be 
expressed as corrected HU values or Tb may be represented in tissue density units such as 
grams per cubic centimeter (g/cc). HU' b is the previously determined or known CT 
number of representative blood, i.e., the known blood density or HU values of blood 
determined in vitro under similar scanning conditions. Blood samples from 
representative patients are placed in a container and scanned inside a tissue-like phantom, 
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while maintaining temperatures representative of the human body, and scanning 
conditions representative of the clinical scans. HU'b is the measured reference tissue 
density which in this example, is the mode of the histogram of the heart and blood. S is 
the slope of the regression equation measured from the calibration phantom samples. HU 
is the scanner CT number or may be a CT number measure of the target tissue to be 
analyzed. Note that the target tissue may be coronary calcifications, an organ, an 
unknown mass, the edge of a target region or organ. Alternatively, the hybrid calibration 
can be applied to the entire image, thus creating a new calibrated image where all image 
voxels are calibrated to the hybrid reference in a step 123. 

[0047] After the image is calibrated, the target tissue or detail to be analyzed is next 
located in a step 124. The target tissue may be located by search algorithms which make 
use of location, density and shape of the targets. In one embodiment, the target tissue is 
located by using the known phantom location to set a search ROI in which a more 
detailed search occurs. Lung calcifications can be found by prior art methods, such as 
described in U.S. Patent No. 4,922,915, cited above. Using the location of the vertebral 
body, the abdominal aorta can be located within an ROI which includes the aorta. The 
measurement ROI may then be placed around the blood of the aorta and the calcifications 
in the aorta. 

[0048] CT images often have streak artifacts, particularly in fast CT scans of the 
heart, which makes it difficult to detect small details such as calcifications. A low pass 
filter can be applied to the image in a step 125 to reduce the streak artifacts without 
affecting the target details, such as the small calcifications in this example. 
[0049] When the target tissue is located and the image calibrated, a pre-selected and 
quantitative threshold is applied or a quantitative threshold is defined as the edge of a 
region in a step 126. The edge of a target or threshold value is located in a step 127 to 
enable reproducible quantitative measurements of dimensions, volumes or mass in a step 
128. Fig. 4 A shows a depiction of an image detail, for example, a round blood vessel. 
When the detail is imaged in, for example, a CT scanner, the final image may have 12 or 
more bits of gray scale values. All of the gray scale values can not be displayed at one 
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time on a monitor, and further, the eye can not see this many gray levels. The image is 
therefore displayed with windows (number of gray levels) and levels (the central gray 
level of the window). A variety of windows and levels are possible in a relatively 
continuous process. Fig. 4 shows a depiction of an image detail and the voxel intensities 
across the detail. A plot across the image would produce the scan line of Fig. 4, showing 
the scan profile with edges. If the image is displayed on an electronic monitor or filmed, 
the perceived location of the edges to the eye will vary with window and level. 
[0050] Fig. 4 also shows how the apparent location of the detail edges changes with 
display level. When the display level is changed from Li to L2, the apparent or measured 
width (distance) increases from Di to D2. As a result, if the level is changed, the object 
appears to get larger and will be perceived as larger by the observer. If physical 
measurements are made of the size or boundaries of the details, they will change, also 
based on the display settings. If the edge is calibrated as disclosed herein, and a 
quantitative edge defined, the detail can be reproducibly measured and displayed. 
[0051] With a defined quantitative threshold or a defined edge parameter or both, the 
boundaries of a target tissue or the number of voxels or both that meet the threshold 
criteria can be determined accurately. The number of voxels satisfying the threshold or 
the extent of the target details can be summed for neighboring CT slices to create a 
calcium score which has 3D definition under the quantitative threshold. The calibrated 
image and a defined threshold enables reproducible detection without operator location of 
details, such as calcifications within the heart or the aorta. Multiple details within a CT 
slice or the volume can be analyzed by repeating the steps 110-128, as illustrated by a 
step 129. 

[0052] Calcifications may be detected and quantified in other parts of the body, for 
example, the abdominal aorta. Figure 7 shows a cross sectional image of the abdomen 
with the reference calibration phantom 5 in place. The abdominal aorta shows a calcified 
plaque 3. The kidneys 10 and the abdominal muscle 12 are also seen. Bowel gas 3 is 
also shown at essentially air density. The subcutaneous fat 1 1 is relatively homogeneous 
and is distributed around the circumference of the body. The subcutaneous fat provides a 
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sharp contrast with the surrounding air. The muscle wall can be segmented out from the 
fat to define a 360° region of purely subcutaneous fat. 

[0053] The results of the foregoing method are reported as a parameter, such as mass, 
volume, or a calibrated calcium score in the case of calcium analysis in the coronary 
arteries or aorta. 

[0054] The foregoing calibration procedure can be used with lesser performance on 
images without the calibration phantom, or with an alternative phantom made of 
alternative materials. A variety of plastics simulating muscle or fat may be used as the 
phantom samples. The use of hybrid calibration is therefore not limited to a particular 
reference phantom or a particular way of being used, i.e., simultaneously scanned versus 
independent scanning. 

[0055] The calibration method of the current disclosure may also be used with 
subcutaneous fat or internal air as a calibration reference for the hybrid method. One 
skilled in the art can see that such alternative methods fall within the slope of the 
following claims. 

[0056] This invention may be embodied in other specific forms without departing from 
the essential characteristics as described herein. The embodiments described above are to be 
considered in all respects as illustrative only and not restrictive in any manner. The scope of 
the invention is indicated by the following claims rather than by the foregoing description. 
Any and all changes which come within the meaning and range of equivalency of the claims 
are to be considered within their scope. 
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